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Abstract 

The screened QED collections to the g factor of Li-like ions are evaluated utilizing an effective potential 
approach. The evaluation is performed within the range of nuclear charge numbers Z = 32-92. The results 
obtained serve for improving the theoretical predictions for the g factor of heavy Li-like ions. 
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High-precision measurements of the g factor of H-like carbon and ox 
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gen 



12, 



stimu- 



13h . In particular, 



lated accurate theoretical calculations of this effect 
these investigations provided a new determination of the electron mass to an accuracy which is 
four times better than that of the previously accepted value (see Ref. [ 14] and references therein). 
Moreover, an extension of these studies to higher Z systems could lead to an independent determi- 
nation of the fine structure constant 1 15, 1(3]. Investigation of ions with more than one electron is 



also anticipated in the nearest future. In particular, measurements of the g factor of Li-like calcium 
are currently in progress by the Mainz-GSI collaboration. 

The motivation for studying Li-like ions alongside with H-like ions arises from the substantial 
elimination of the uncertainty due to the nuclear size effects in a specific combination of the cor- 
responding g factor values II 1711 . The most accurate results for different contributions and for total 
values of the g factor of Li-like ions within the range Z = 6 - 92 were presented in our recent work 
ifisl . For all values of Z the uncertainty of the total value was mainly defined by the contributions 
of the interelectronic interaction and the QED screening effect. The latter was evaluated to the 
leading order (a/n)(aZ) 2 /Z and calculations to all orders in aZ are needed for high values of 
Z. In the present paper we perform this evaluation using an effective potential approach. The 
obtained results are combined with other contributions to improve theoretical predictions for the 
g factor of Li-like ions within the range Z = 32-92. Compared to our previous work [ 18], the 
new theoretical values include also recent results for the higher-order two-loop QED corrections 
J12I and the magnetic-loop part of the vacuum polarization Jisl . In addition, the nuclear-size 
correction is recalculated employing the most recent data for nuclear radii from Ref. [19]. 

The relativistic units (h — c — m — 1) and the Heaviside charge unit (a = e 2 / (4.7r), e < 0) are 
used throughout the paper. 

The total value of the ground- state g factor of a Li-like ion is conveniently written as 

g = g D + Atfcnt + As-qed + A# S qed + A# nuc , (1) 

where 

g D = l(l + V /2 + 2 v /r^Z^) = 2 - M£ + . . . ( 2 ) 

is the Dirac value for the point-charge nucleus, Ag int is the interelectronic-interaction correction, 
A#qed is the one-electron QED contribution, A(?sqed is the screened QED correction and Ag nuc 
incorporates the nuclear- size, nuclear-recoil and nuclear-polarization corrections. For the eval- 
uation of the interelectronic-interaction and nuclear corrections we refer to our previous papers 



OH 



and references therein. Radiative corrections for electrons bound in a pure Coulomb field 
were addressed, in particular, in Refs. . The main goal of the present work is to 

evaluate the screened QED corrections within an effective potential approach. 

We consider an effective spherically symmetric potential V e g that partly takes into account the 
interelectronic interaction between the valence 2s electron and the core electrons of the (Is) 2 shell. 
The simplest choice of V c r is the core-Hartree (CH) potential 

f°° 1 

V eS (r) = Kuc(r) + a / dr'—p c {r'), (3) 



r- 



ear potential. The screening 



where p c is the density of the core electrons and V^ uc denotes the nuc 

potential derived from the density-functional theory reads J20L 21II22" 

poo j a / 81 \ 

V e $(r) = V nuc (r) + a dr' — p t (r') — x a ( — — ^(r) . (4) 
Jo r> r \32n 2 J 

Here p t is the total electron density, including the (Is) 2 shell and the 2s electron. The parameter 
x a can be varied from to 1. The cases of values x a = 0, 2/3 and 1 are referred to as the Dirac- 
Hartree (DH), the Kohn-Sham (KS) and the Dirac-Slater (DS) potentials, respectively. To provide 
a proper asymptotic behavior, equation © should be replaced by H23I1 

r 

at large r, where N c is the number of core electrons. The self-consistent potential @ is generated 
by iterations, which continue until the energies of the core and valence states become stable on the 
level of 10~ 9 . The CH potential © does not imply self-consistency and, therefore, it is generated 
after one iteration. When the effective potential V c q and the spectrum of the corresponding Dirac 
equation are generated, the screening correction A(?sqed is calculated as the difference between 
the two values of the QED correction calculated for the potential V c s and for the nuclear potential 

'nuc 5 

A#SQED = A^qed [V cS ] ~ A#qed [Kuc] • (6) 

Below, we describe the evaluation of the one-loop QED correction to the g factor in arbitrary 
binding potential. The QED correction of first order in a appears as the sum of self-energy and 
vacuum-polarization corrections, A^q^ d = Ag SE + Ag V p- The vacuum-polarization term is 
relatively small, and its contribution to the QED screening can be neglected at the present level of 
accuracy. The self-energy correction is given by the sum of irreducible, reducible and vertex parts, 

A#se = A5- irr + Ag red + Ag vei . (7) 



The expression for the irreducible part reads H24fl 



= J_ £ \p (aj (S(gg) - 7°^) \n)(n\[r x g^g) 



Here |a) is the 2s state with the angular momentum projection m a , 5m is the mass counter-term 
and E(e) denotes the unrenormalized self-energy operator defined by 

27ri-oc ^ e - w - e„(l - iO) 

where J (a;, x x , x 2 ) = e 2 a ll oi u D tlv {uj, xi — x 2 ), ct M = (1, a) are the Dirac matrices and is the 
photon propagator. To separate ultraviolet divergent zero-potential Agr^ and one-potential Agr^ 
terms we follow Ref. H25I1 and calculate them in momentum space. The residual part of Ag irr , 
so-called many-potential term Agf; rr ' , is calculated in coordinate space employing the algorithm 
proposed in Ref. [26]. The expression for the reducible part reads 12411 

1 d 

A&ed = («k- S ( £ )L=e l a )( a l[ r X <x]z\a) , (10) 

m n oe a 



while the vertex part is given by [24] 



A, ver = ±±- du Y (^l^)|nia)K|[rxa]> ; ) 

m a 2ir y_ 0O ^ (e a - u - e ni (l - i0))(e a - a; - e„ 2 (l - zO)) 



Both reducible and vertex parts are ultraviolet divergent, whereas the sum Ag vr = Ag ver + Ag red 




is finite. Following Refs. (5|,lUJ], we separate out the zero-potential term Agi? and evaluate it in 
momentum space. The remaining many-potential term Aglr + ' ) is calculated in coordinate space as 
the point-by-point difference between the contributions (ITU with bound and free propagators in 
the self-energy loop. 

Below we discuss some details, concerning coordinate- space calculations of A^ r 2 r + ^ and 
Agil + \ Angular integration and summation over intermediate angular momentum projections 
is carried out in the standard way. The many-potential terms involve infinite summation over the 
angular quantum number k. The summation over the complete spectrum of the Dirac equation at 
fixed k is performed using the dual-kinetic-balance approach [27]. The summation over k is termi- 
nated at \k\ = 10-20 and the rest of the sum is estimated by the least-square inverse-polynomial 
fitting. It was observed that more stable values are obtained, when subtraction © is performed 
prior to the fitting procedure. 
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The results obtained for the screened QED correction are presented in TableU The uncertainty 
of the numerical evaluation is defined by the many -potential terms. The uncertainty due to incom- 
pleteness of the effective potential approach in the description of the interelectronic interaction is, 
however, much larger. It was estimated as the mean deviation of the results obtained by means 
of the different potentials: core-Hartree, Kohn-Sham, Dirac-Hartree and Dirac-Slater. Since for 
middle-Z ions the accuracy of our previous results 111 8ll turns out to be better than that of the 
present ones, in the final compilation (see below) for Z < 30 we present the values of A^sqed 
from 1 181. We mention that these values include terms of higher order in 1/Z calculated by Yan 

n 

[28J. For Z > 30 we take the present data, obtained with the Kohn-Sham potential. While for 
middle Z there is a good agreement between the present and the previous results, a significant 
discrepancy is found for high values of Z. This is presumably due to the fact that the method of 
Ref. II 1 811 - based on the non-relativistic treatment of the anomalous magnetic moment, is rather 
insensitive to the bound- state QED effects. The accuracy of the present evaluation of the screened 
QED correction is limited by the possibility to account for the interelectronic-interaction effects in 
terms of the local potential V c s. A rigorous evaluation of the 1/Z contribution to A^sqed should 
be the next step for improvement of the accuracy of the many-electron QED correction. 

In Table UH the individual contributions to the g factor of the ground state of Li-like ions are 
presented. The changes made compared with the previous compilation lisl concern several terms: 
the one-electron QED correction, the screened QED correction and the finite-nuclear-size correc- 
tion. The evaluation of A^sqed is already described above. The finite-nuclear-size correction is 
recalculated with the most recent data for the nuclear radii ifioll . The one-electron QED correction 
of first order in a is updated with the recent evaluation of the magnetic-loop vacuum-polarization 
term 111311 . This reduces the uncertainty of the A^qed term for Z > 30. The one-electron QED 
correction of second order in a is improved for Z < 30 employing the analytical formula for the 
(a/7r) 2 (aZ) 4 term recently derived in Ref. I12I1 . 

In summary, we have presented the evaluation of the screened QED correction to the g factor 
of Li-like ions within the effective potential approach. These results improve the accuracy of the 
theoretical predictions for the g factor within the range Z = 32-92, where stringent tests of the 
bound-state QED effects are expected. More elaborated treatment of the interelectronic-interaction 
and screened QED corrections will be the subject of our subsequent investigations. 
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TABLE II: Individual contributions to the ground-state g factor of Li-like ions. 
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